arXiv:1504.04978v2 [cond-mat.mes-hall] 8 May 2015 


Longitudinal conductivity of massless fermions with tilted Dirac cone in magnetic field 


Igor Proskurin/’^’EI Masao Ogata/ and Yoshikazu Suzumura^ 

^Department of Physics, The University of Tokyo, Bunkyo-ku, Tokyo 113-0033, Japan 
^Institute of Natural Sciences, Ural Federal University, Ekaterinburg 620083, Russia 
^Department of Physics, Nagoya University, Nagoya 464-8602, Japan 
(Dated: May 11, 2015) 

We investigate a longitudinal conductivity of a two-dimensional relativistic electron gas with a 
tilted Dirac cone in magnetic field. It is demonstrated that the conductivity behaves differently in 
the directions parallel and perpendicular to the tilting of the cone. At high magnetic fields, the 
conductivity at non-zero Landau levels in the direction perpendicular to the tilting modifies non- 
trivially, in contrast to the parallel case. At zero temperature, the crossover of the conductivity at 
the Dirac point from high to low magnetic field is studied numerically. It is found that the tilting 
produces anisotropy of the conductivity which changes with the magnetic field which is different 
from the anisotropy coming from the Fermi velocity. We also discuss the conductivity at finite 
temperatures and finite magnetic fields which can be directly compared with the experiments in a- 
(BEDT-TTF) 2 l 3 organic conductor. We find that the tilting does not affect so much the magnetic- 
field dependence of the conductivity except for the prefactor. We discuss the interpretation of recent 
experimental data and make some proposals to detect the effect of the tilting in future experiments. 

PACS numbers: 72.10.-d,75.47.-m 


I. INTRODUCTION 

During the last decade, Dirac fermions have gained 
much attention in condensed matter physics being of a 
great interest from both fundamental and applied points 
of view.d. The most prominent material where existence 
of massless Dirac particles was unambiguously shown 
is graphenei^ Later Dirac electrons were theoretically 
predicted and experimentally observed in many other 
materials including surface states of two- and three- 
dimensional topological insulators^ a-(BEDT-TTF) 2 l 3 
organic conductor^^^S and BaFe 2 As 2 iron-pnictidei^ The 
presence of relativistic Dirac quasiparticles in these ma¬ 
terials gives rise to such interesting physical phenomena 
as unconventional quantum Hall effect;^ crossover from 
the positive to negative interlayer magnetoresistance) ^ ^ 
and giant Nernst-Ettingshausen effectji2rJ4 

Organic conductor q:-(BEDT-TTF) 2 I 3 is a layered ma¬ 
terial which consists of conducting layers of BEDT-TTF 
molecules separated by insulating layers of I;^ anions. 
Under the hydrostatic pressure about 1.5 GPa, it un¬ 
dergoes a transition to the zero-gap state with a typ¬ 
ical ratio of the in-plane to the interlayer conductiv¬ 
ity ~ 10 ^ which makes this compound a quasi-two- 
dimensional conductor Unusual transport properties^ 
in a-(BEDT-TTF) 2 I 3 were explained by finding the ex¬ 
istence of the Dirac fermions in the analytical band 
structure calculation,- and numerical onei^ Kobayashi 
et alX showed that Dirac quasiparticles should be de¬ 
scribed by an anisotropic tilted Weyl equation. Fur¬ 
ther experimental and theoretical studies provided con¬ 
firmation of the massless Dirac nature of the carriers 
in this compoundiIn contrast with graphene, 
where the carrier concentration is easily tuned by apply¬ 
ing the gate voltage, the carrier number in a-(BEDT- 
TTF) 2 l 3 is always fixed near the charge neutrality point 


which makes difficult direct observation of the Landau 
levels. Recently, multilayer quantum Hall effect and 
Subnikov-de Haas oscillations were observed in this com¬ 
pound using the method of contact electrification from 
the substrate 

A notable feature of the massless fermions in a- 
(BEDT-TTF) 2 l 3 is considerable tilting of the Dirac cone 
in the energy-momentum space which breaks both the 
Lorentz invariance and the chiral symmetry. The prop¬ 
erties of systems with a tilted Dirac spectrum and their 
realizations in real materials were studied by several 
authorsJ^ 4 ^ It was found that in magnetic field the 
effect of tilting on the energy spectrum is similar to 
the effect of an electric field on relativistic Landau lev¬ 
els in graphene.i^Si^i In the latter case, using the ar¬ 
gument based on the Lorentz covariance of the Dirac 
Hamiltonian, it was shown that the electric field squeezes 
the separation between the Landau levels and shifts 
the positions of the wave functions in the real spacej^ 
Kawarabayashi et al^ found that even if the tilting is 
introduced, there exists the generalized chiral symmetry 
which provides existence and stability of the zero-mode 
which governs the transport properties near the charge 
neutrality point in high magnetic fields 

A general problem of the transport of massless Dirac 
fermions was studied previously in different contexts. In 
relation to graphene, Shon and Ando2^ calculated the 
conductivity in the self-consistent Born approximation 
(SCBA), and Peres et al^ considered the effect of lat¬ 
tice defects and electron-electron interactions. Sharapov 
et al^ made a fully analytical calculation of the zero- 
temperature conductivity for gapped Dirac excitations 
in the context of d-wave superconductors, and the prob¬ 
lem of Shubnikov-de Haas oscillations was studied in de¬ 
tail by Gusynin and Sharapov With respect to organic 
conductors, the effect of electron-electron interactions on 
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the in-plane conductivity was considered by Morinari and 
Tohyamaj^ The effect of tilting on the in-plane conduc¬ 
tivity in zero magnetic field was studied previously by the 
authorsi^SiSi In weak magnetic fields the effect of tilting 
on the Hall conductivity has been studied previously by 
Kobayashi et al^ 

The purpose of this paper is to answer several ques¬ 
tions. The first question is how the tilting of the Dirac 
cone affects the transport properties of two-dimensional 
Dirac fermions in magnetic field. The second question 
is how to provide an experimental confirmation of the 
tilting. According to the knowledge of the authors, at 
the present time, there is no experimental evidence of 
the Dirac cone tilting in a-(BEDT-TTF) 2 l 3 . In this pa¬ 
per, we make a proposal how it can be confirmed in fu¬ 
ture experiments by detecting the anisotropy in the con¬ 
ductivity. The third question is related to the scenario 
of the magnetotransport in a-(BEDT-TTF) 2 l 3 which is 
motivated by recent experimentsi^^*^ We note that at 
the present time the magneto-conductivity in a-(BEDT- 
TTF) 2 l 3 is not completely understood. First-principles 
band-structure calculations show that, together with the 
Dirac cone, a-(BEDT-TTF) 2 l 3 also has a heavy-hole 
band near the Fermi level^ which might participate in 
the conducting properties. Based on this argument, the 
magneto-conductivity of a-(BEDT-TTF) 2 l 3 in Ref. [H 
was explained in the framework of quasi-classical two- 
carrier model where the first carriers are massless Dirac 
electrons, and the second are massive holes. However, 
taking into account recent experimental results^, we pro¬ 
pose an alternative explanation which is based only on 
the Dirac type of carriers. We discuss a criterion which 
can make a clear distinction between these two models. 

The paper is organized as follows. In Sec. m we de¬ 
velop the formalism. In Sec. uni we consider the ef¬ 
fect of the Dirac cone tilting on the zero-temperature 
conductivity. In the limit of high magnetic fields, we 
study how the tilting modifies the conductivity when the 
chemical potential is situated at the nth Landau level. 
It is demonstrated that the conductivity in this case 
is determined by the wave function of the nth Landau 
level. Next, we study how the conductivity at /i = 0 
depends on magnetic field and impurity concentration. 
From the previous studies^^i^i we know that massless 
Dirac fermions at /r = 0 possess universal conductivity 
which is magnetic field and impurity independent. We 
show that in the presence of the tilting this behavior, 
in general, changes, however, in one particular direction, 
the conductivity remains independent of impurities and 
magnetic field. We also explain that the tilting leads 
to the anisotropy in the conductivity which changes in 
magnetic field, which makes it different from another 
type of the anisotropy coming from the Fermi velocity. 
This fact can be useful for experimental confirmation of 
the tilting in the two-dimensional Dirac systems. Sec¬ 
tion E] is devoted to the conductivity at finite temper¬ 
atures. First, we demonstrate that a combination of a 
crossover from weak to high magnetic fields, and the 


Zeeman splitting of n = 0 Landau level result in char¬ 
acteristic two-step decrease of the conductivity in the 
magnetic field detected experimentally)^!^ Next, based 
on the self-consistent Born approximation (SCBA) argu¬ 
ments, we show that the magnetic-field dependence of 
the Landau level broadening produces minimum in the 
magneto-conductivity which may be crucial for under¬ 
standing recent experiments)^ We also find that the tilt¬ 
ing does not affect so much the magnetic-field dependence 
of the conductivity except for the prefactor. Section IVl is 
reserved for the summary. 


II. FORMULATION 

A. Model Hamiltonian 

The Hamiltonian for the two-dimensional relativistic 
electron gas inside the layer of BEDT-TTF molecules in 
magnetic field (for a given valley) can be written in the 
form of the generalized Weyl Hamiltonian^iiS 


H = Vf {ll^yCro + ^xCTx + IlyCry) , (1) 


where vp is the Fermi velocity, H^ = —ihVi + eAi denotes 
a canonical momentum, Ui is the Pauli matrix (i = x, y), 
with (To being the unitary matrix, Ai is a magnetic vector 
potential, and —e is an electron charge. Here, we consider 
the case of a tilted Dirac cone with the isotropic Fermi 
velocity due to the fact that the angular dependence of 
the Fermi velocity is small in Q!-(BEDT-TTF) 2 l 3 )^ We 
comment on the effect of the anisotropy of the Fermi 
velocity at the end of Sec. Hill Without loss of generality, 
we imply that the Dirac cone is tilted in the y-direction 
with 0 < y < 1 being a degree of the tilting. 

The eigenvalue problem for the Hamiltonian in Eq. o 
in the Landau gauge A = (0, Bx, 0) can be solved by al¬ 
gebraic methods which gives the following spectru m^^i^^ 


En = sgn{n)^j2eBhvjp\^\n\, 
and eigenfunctions 


'^knix,y) = 


giky 

VeL 




X — X„ -I- i'^k 


( 2 ) 

(3) 


where 


- isgn(n)x^+^fe|n|-i(x) 

\/2(2 — i5„o)(l + A) 


(4) 


n = 0, ±1, ±2,..., A = £ = \Jh/eB is the mag¬ 

netic length, k and L are the wave number and the size of 
the system in the y-direction, where peri odic bo undary 
conditions are implied, x„ = — sgn{n)r]£\j2\n\/\, = 

(1 -I- A, —iri)^ = (*^) 1 + ^)'^] denotes the eigenvec¬ 
tor of a generalized chiral operator— 7 = (uz — irjCJx) /A 










3 


corresponding to +1 (—1) eigenvalue, 

’ 

(5) 

and Hn{x) denotes the nth Hermite polynomial; sgn(x) 
is defined as —1 for x < 0, 0 for a; = 0, and +1 for x > 0. 
Each Landau level is multiply degenerated with respect 
to k with the degeneracy factor E/(27r£^) where E is a 
two-dimensional area of the system. 

In what follows, we consider the Hamiltonian in Eq. m 
for a given valley index and spin projection. The effect 
of the Zeeman interaction between the electron spin and 
the magnetic field can be taken into account by including 
the term —gT^fJ-sB/2 where Tj, is the Pauli matrix of the 
real spin, /tb is the Bohr magneton, and g denotes the 
^-factor. For the sake of brevity, we will omit the Zeeman 
term in most of the formulas following and restore it when 
it is necessary. 


B. Longitudinal conductivity 


At zero temperature, the longitudinal part of the con¬ 
ductivity tensor for an electron gas in the presence of 
static disorder can be calculated using the Kubo-Bastin- 
Stfeda formul a^^i^'^ 

aii{0,g) = — (Tr [xj ImG'(/r)x 2 ImG'(^)]), (6) 

TT 

where ^ is a chemical potential, the velocity operator is 
defined as Vi = (i/h) [H^n], ImG = (i/2) — G^”^], 

and G(^^(e) denotes {e — 'H± i6)~^. Here, % = Lf-I-Idmp 
is a total Hamiltonian which includes an impurity poten¬ 
tial Idmp (specified later). The angle brackets stand for 
the average over the impurity positions. As was shown 


previously, in the case of short-ranged scatterers the cur¬ 
rent vertex corrections in Eq. m vanish and here¬ 
after we imply that the effect of the impurities is included 
into the Green function self-energy. The conductivity at 
the finite temperature can be restored from the expres¬ 
sion 


au{T,n) = -f de^^^a{0,e) (7) 

where /(e) = {1-I-exp [(e —/r)//cBr]}~^ is the Fermi- 
Dirac distribution function with the temperature T. 

Calculating the trace in Eq. dH]) with respect to the 
Landau level eigenfunctions we obtain 

„2t 2 °° 

\{n\vi\m)f An{n)Am{fi), (8) 

m,n— — oo 


where i = x,y, 


and 


Vx = cTx and Vy = ay + g, 


An (/^) 


__ 

ifM-En-Rig)f + T^fiy 


(9) 

( 10 ) 


We imply that the electron self-energy E(^) = R{fi) + 
is free of the Landau level index dependence. Here, 
an denotes the conductivity per valley and per one spin 
projection. The Zeeman interaction in Eq. m can be 
easily restored by the substitution aii{0,fi) —>■ aii{0, g + 
-I- aii{Q,g- ^gsB). 

The matrix elements of the velocity operators can be 
calculated explicitly 

(m\Vx\n) = -iXAn-mP^igAn,m), (11) 

{m\Vy\n) = An, {g An,m), (12) 

where n ^ m, An,m = sgn(n)y^- sgn(m)y^|mi. 


P:^{j 1 An,n.) = 


^{2-5no){‘i-5mo) 




X I '’"'(??^A2_„)-sgn(n)sgn(m)^/T^Lj()||J)"'(??2A2^^) J , ( 13 ) 


and L'^(x) denotes the generalized Laguerre polynomial. 
The matrix elements of velocity operators with n = m 
equal to zero. The details of the derivation are given in 
Appendix [A] 

For simplicity of the subsequent analysis, we introduce 
dimensionless conductivities in the direction perpendicu¬ 
lar a± = ax:x/ and parallel a\\ = ayyj to the 


tilting where h = 27rft. The dimensionless conductivities 
can be written as functions of dimensionless parameters 

= ^n,-m {P^iV^n,m)f gm(x)g„(x),(14) 

n,m 

^11 (^) = ^ ^ ^n.m \Bni {V^n,m)] 9m(,x)gn{x), (15) 

n,m 











4 


where 


gn{x) 


_7_ 

{x - sgn(n)y ^)2 + 7 ^ 


(16) 


X = {g- R{fJ,))/{hw*), 7 = and fuv* = 

is a cyclotron energy huic = \/2ehvpB renor¬ 
malized by the tilting. 

In the r] — 0 limit the velocity operators in Eqs. m 
and m have matrix elements only for Landau levels 
with |n| = \m\ ± 1 , and the expressions for the conduc¬ 
tivities in Eqs. m and m are considerably simplified. 
The conductivity becomes isotropic = ux = where 


_ nuc ^ ^_ 7 _ 

° 4 “ a-v/rr+T)2 -|- 72 

X - -- 7 . (17) 

{x — a'y/n)^ + 72 

The summation over the Landau levels in this case can 
be done analytically^^ which gives 


ao{x) = 1 - 


x'^ + 7 ^ 

1 -I- 16a;^72 


x^ + 4x^7^ -I- 8x^"f^{x^ + 7 ^) 


(a:^ -I- 72)2 


-4x7lm'!/i ^(7 -I- ixY 


(18) 


The details are given in Appendix [BJ In the case when 
i?(/i) = 0 and r(p,) is a constant (constant damping 
approximation), the conductivity at low magnetic field 
can be approximated using the expansion ip{z) = log z — 
l/(2z) — l/(12z^) -I- ... which yields 


^XX (Oj 1 ^) 



32r4 Hr/ ’ 


(19) 


where <i)o(a:) and $ 2 ( 2 ;) are defined by 


4>o(a;) 


1 

2 




( 20 ) 


and 


$ 2 ( 2 :) 


1 — 

1 + x"^ 


1 -b a;^ 

X 


tan ^ X 


8a:^(l — x^) 
^ 3(l + a;2)3 


( 21 ) 


The expression for $0 is well knowni2^i2i In the |a;| ^ 1 
limit, we obtain $0 ~ ■n\x\/2, $2 ~ 7 r/( 2 |a;|), while $0 = 
1 -b a;^/3, $2 = 128a;^/15 for a; <C 1. 


C. Self-consistent Born approximation 

In order to take into account the effect of scattering on 
the impurities in a self-consistent manner, we consider 
the total Hamiltonian in the form of "H = iL -b Idmp where 

Mmp = 5 (r-i2,) (22) 

3 


NsU^ 



Figure 1. Diagrammatic representation of Eq. (|23l). Each 
dashed line on the right-hand side stands for uS(r — Rj) and 
the star means the average over Rj and summation over j. 


describes the interaction between the conduction elec¬ 
trons and randomly distributed point-like impurities at 
the positions Rj with the scattering potential u. Here, we 
imply that the impurity potential does not mix the val¬ 
leys, which corresponds to the case of short-ranged scat- 
terers in Ref. [ 2 ^. The self-energy E(e) in SCBA can be 
calculated using the method developed by Bastin et al^ 
and Shon and Ando2^ which leads to the self-consistent 
equation S(e) = Nsu‘^V~^ Ti G{e) illustrated in Fig. [T] 
The explicit form for this equation reads 


S(e) 


27r€2 ^ e_£;^_E(e)' 

n— — oo ^ ' 


(23) 


where A^s is the impurity concentration. Here, we ignore 
the effect of the Zeeman interaction for simplicity which 
is discussed below. 

In zero magnetic field, we can replace the summation 
over the Landau levels in Eq. (E21) by integration. In 
the Boltzmann limit, where we neglect the real part of 
S(e) and suppose that r(e) = ImS(e) is negligible, we 
obtain^^ 

r(e) = ^|e|, (24) 

where k = NgV?'/{27rvph?) is a dimensionless parameter. 
The only difference of Eq. (El from the case of 77 = 0 is 
the appearance of the factor which comes from the 
angular dependence of the density of states at the Fermi 

leveliSi^ 

H 1 1 _ 1 _) 

Jq 27r (1 — ?7sin(/))2 _ ^^2^3/2 ^3 

In the opposite limit, when e —>■ 0 first, we obtain 
r(0) = Ac exp [—A^/( 2 >^)] where we introduced the cut¬ 
off energy to regularize the logarithmic divergence. 
The only difference of this result from 77 = 0 case studied 
in Ref. is the factor A^. 

In the quantizing magnetic field, when e E^, we can 
keep in Eq. (1231) only the contribution with n = N. In 
this case, r(e) is approximated as 2 ^ 


r(e) 





2 (e-Ajv)2 




K 


(26) 
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Figure 2. (Color online) ImE(e) at e = 0 as a function of hcuc 
for several k with Ec = 1 and 77 = 0. 


energy at e = 0. In order to obtain Fig. O we solved 
Eq. (l!I(Tll numerically with fixed Ec = 1. At Huic = 0 the 
value of Im E( 0 ) for small becomes exponentially small 
in agreement with the expression exp [—A^/(2>r)]. For 
large enough cyclotron energy, the dependence becomes 
linear which corresponds to r(0) ^ In the follow¬ 

ing, we will denote 

r(B) = ImE(O) - (33) 

obtained in SCBA in order to distinguish it from the 
constant F. 


III. CONDUCTIVITY AT ZERO 
TEMPERATURE 


In finite magnetic fields, in order to have a convenient 
analytic expression for Eq. (1231) . we rewrite it in the fol¬ 
lowing form: 

E(e) = r du (27) 

J_oo e-a;-E(e) 

where 

1 

^oH = ^ E (28) 

^ n—— 00 

is the density of states in the absence of the impuri¬ 
ties. Equation dia contains a divergence which can be 
treated analytically using the convenient representation 
for Do{oj) obtained in Ref. IH: 


In this section, we analyze the effect of the tilting on 
the zero temperature conductivity in high magnetic fields 
where the Landau quantization of the energy spectrum 
plays the principal role. At fist, we analyze the chemical 
potential dependence of the conductivity using Eqs. (O 
and m- After that, we give a clear physical explana¬ 
tion of the conductivity at the Nth Landau level and 
discuss its relation for the quasi-classical picture. Sec¬ 
ond, we numerically study how the conductivity at the 
Dirac point depends on magnetic field and impurity con¬ 
centration. The results are given for two models of the 
impurity scattering: constant broadening approximation 
and SCBA. 


A. High magnetic fields 


Do{uj) = 


1 

H- — tan 

TT 






cot 


{hiui 


■ (29) 


Now the integration in Eq. (EH) can be performed explic¬ 
itly which finally gives the self-consistent equation 


where 

$(a;;ec) 


E 


= X [log (e 


A3 \ hiu* ^ 


1 

~ 



TT cot(7ra:^) 


(30) 


(31) 


The effect of the Zeeman splitting can be introduced in 
Eq. (1501) by changing the right-hand side to 





— E — 



(32) 


We start with analyzing the conductivity at T = 0 
in the high magnetic field when Landau levels are well 
separated ihiUc ^ r(/i)). In this case, when the chemical 
potential is close to Ejy we can keep only the terms with 
n = N in the summation in Eqs. da and da which 
gives 




_ ^{1} _ I"w(m) _ 


(34) 


where I = (T,||). The coefficients which are 
at /i = Ajv, are defined by 

= Wx E [P^^ivAN,m)y , (35) 

[PN^ivAN,mf. (36) 

m^N 

For the constant broadening approximation the conduc¬ 
tivity in the vicinity of the TVth Landau level is given 

by 


The e dependence of ImE(e) was studied in Ref. [ 2 ^. 
In Fig. [21 we show how Im E(e) depends on the cyclotron 




( 37 ) 
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while if we use Eq. (uni) for r^v in SCBA, the expression 
for the conductivity becomes 




XT [huJcY 


(38) 


Note that ct® does not depend on T{fi), and the differ¬ 
ence between constant F and SCBA appears only in the 
/i-dependence of 

The ry = 0 case was studied by Shon and Ando^^ who 
showed that cry = + 2|iV|). For finite rj, using 

the numerical summation over Landau levels in Eqs. dSSD 
and (ESI), we have found similar results 

^ {Sm + m\ - Nif) , (39) 

= ^^{5m^ + 2\N\). (40) 

Note that for TV = 0, these equations can be verified 
explicitly using Eq. o, 




lyii) 


[rj’^nY 




-e-" " = 1. (41) 


while for \N\ > 1 the expansions (IB7I) and (IBllI) allow 
to confirm these equations up the order of rj^. 

From Eqs. (El and (HiH) we find out that both crj_ and 
cr|| are modified by the tilting through the factors 
and A, respectively, which can be accounted for the geo¬ 
metrical modification of the quasi-classical electron orbit 
(discussed below). However, crj_ has an additional mod¬ 
ification by the tilting for N ^ 0 a.s demonstrated in 
Fig. [21 where the normalized conductivities cr|| = A^^cry 
and crj_ = Acrj^ are shown as functions of /i calculated nu¬ 
merically in SCBA using the summation over 50 Landau 
levels in Eqs. m and m- 

The reason why cr|| at /i = remains the same as 
in the 77 = 0 case, while a± is reduced by Nrj'^/X can be 
understood as follows. From the Eq. the conductivity 
at ^ = En can be expressed as 


a 


(0 _ 
N ~ 


e^h ^ {N\vi\m) (m\vi\N) 
{En - E^f 


(42) 


where the prime indicates that the m = N term is omit¬ 
ted. Using the relation vi — {i/K)[E[,ri]^ we can write 
the off-diagonal matrix elements of the velocity operator 
as 



p/fiw* 

Figure 3. (Color online) Chemical potential dependence of 
the normalized conductivities cth = A~^cr|| and = Actx in 
SCBA with 77 = 0.5, x = 0.005, and Ec/hu)l = 100. The 
dotted-dashed (dotted) horizontal line is a guide for eyes at 
the value 2N (2A — Nrf). 


at fj, = Eff in the quantum limit is determined only by 
the Nth Landau level wave function and proportional to 
the quantum-mechanical average of (Ar/)^ in the 
state. 

From the semiclassical point of view, the trajectories 
of a Dirac electron in the momentum space with the 
tiling 0 < 77 < 1 are displaced ellipses which can be 
parametrized as 


kx = bN sin (j), (45) 

ky = xn£~‘^ + ttN cos (j), (46) 

with the semimajor (semiminor) axis equal to a n = 
i-^^y2\N\/\ {bN = and 0 < (/> < 27r. The 

area of each ellipse is independent of rj and given by 
Sn = 27r|7V|t'“^. Note that the Berry phase contribu¬ 
tion to the semiclassical quantization rule = tt, as 
in the case without tiltingJ^ Each ellipse has its focus 
at the axes origin and its center displaced by in 

the fcy-direction, as shown in Fig. The eccentricity of 
each ellipse is rj. The trajectory in the real space can be 
obtained by rotation of the orbit in the momentum space 
by 7r/2 and rescaling it with the factor If we make a 
semiclassical calculation of {{Ari) )n, we obtain 


(iV|7;/|m) = ^ {En - Em) {N\ri\m) 
n 

which gives 

p) _ {{Anf)N 
^ TT^h e 


(43) 


(44) 


where {{AriY)]sf = (N\rf\N') — (7V|n|iV)^. Even in the 
presence of the tilting, the zero-temperature conductivity 


{{Axr)N = fj g (X - x^f = (47) 

and ((Ay)^)Ar = A|A|. The appearance of the factor A“^ 
(A) in the ctx (ch) is therefore accounted for the elliptic- 
ity of the electron orbit. However, in the semiclassical 
picture the displacement of the orbit has no effect on the 
conductivity and behavior of ax and cr|| at fi = En is the 
same as in the case of 77 = 0 . 
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Figure 4. (Color online) Semiclassical trajectories of a mass¬ 
less Dirac electron with the tilted cone in the k space with 
77 = 0.6 for n = —5,... ,5. The units of k^ and ky are £~^. 
The blue (purple) ellipses correspond to the upper (lower) 
band. The dots show centers of the ellipses. Each ellipse has 
its focus at the origin and its center displaced by Xn£~^ in the 
7/-direction. The eccentricity of each ellipse is t/. 


The quantum-mechanical average of the and ri 
{I = X, y) operators calculated with respect to the eigen¬ 
functions in Eq. (ED gives the following results: 


{N\x\N) = ^XN, (48) 

{N\x^\N) = 2xI + ^{5no + 2\N\), (49) 

(iV|y2|iV) = —(5jvo + 2|iV|), (50) 


and {N\y\N) = 0 which justifies Eqs. (IMl) and (l40l) ob¬ 
tained by numerical calculations. These results are dif¬ 
ferent from those obtained in the semiclassical picture, 
because (A^|a;|iV) ^ xn- Note that according to our 
gauge choice, y is an unbounded operator and a regu¬ 
larization procedure is required. The details are given in 
Appendix [Cl 


B. Magnetic-field dependence 

It is well known that in the rj = 0 case, the zero- 
temperature conductivity at /r = 0 is independent of 
magnetic field and impurity concentration and given by 
the universal value The natural question is 

how this behavior is modified for 77 > 0. At zero mag¬ 
netic field, we already know the answer. In our previous 
work, we have found that the conductivities at /i = 0 are 
independent of the impurity broadening and defined by 


the following expressions)^ 


1 

r, - 2 ’ 

vi-r 

(51) 

sin“^ 77 
"^11= ^ ■ 

(52) 


If we compare these results with the case of strong mag¬ 
netic field given by Eqs. (1571) and (1551) with iV = 0, we 


would find that (Tj_ remains the same in both limiting 
cases huj* —>■ 00 and fiw* —>■ 0 , while a\\ increases from 
— 77 ^ < 1 to sin“^ 77/77 > 1 as hu* is reduced from 00 
to zero. 


Eirst, we analyze the behavior of the conductivity at 
low magnetic fields in the constant damping approxima¬ 
tion. In order to study the behavior of the conductivi¬ 
ties in the region of moderate magnetic fields, we have 
made a numerical summation in Eqs. (HI and m over 
N = 150 Landau levels. Eigure [5] (a) shows normalized 
conductivities d-± = Xa± and = A“^cr|| at /r = 0 for 
moderate tilting 77 = 0.4 and 0.5. Note that at 77 = 0 
the only parameter is E /fyxi* which is proportional to 
. Within the numerical accuracy we found that (t_l 
remains constant with increasing T/ hio* and does not de¬ 
pend on rj, while (T|| increases with increasing T or 77, but 
remains smaller than the corresponding limiting value 
sin“^ ri/{\ri) shown by dashed lines. Note that for large 
values of 77, numerical calculations become complicated 
especially due to the inter-band matrix elements of veloc¬ 
ity operators with large |m| and |n| — | 7 T 7 | [see Eq. USD]- 


Second, we show the normalized conductivities in 
SCBA in Fig. |5][b) where Eq. (1551) has been used. For 
the self-consistent calculation we kept Ec/hoj* fixed and 
varied parameter >c which is proportional to the total 
number of impurities. In SCBA, the behavior of the con¬ 
ductivities in the crossover region remains qualitatively 
the same as discussed above. As >c increases, within 
the numerical accuracy cf^ remains constant, while 
increases, but remains below the limiting value at zero 
field (which is the same in both constant broadening and 
SCBA cases). 

At the end of this section, we comment about the 
difference between the anisotropy of the conductivity 
given by the tilting and the anisotropy due to the dif¬ 
ference of Fermi velocities Vx ^ Vy in the x and y- 
directions. In the latter ca se, the energy spectrum iaiS 
En = sgn{n)vp^2eBh\n\, where Vp = ^VxVy. The 
quasi-classical orbits are ellipses centred at the axes ori¬ 
gin with the eccentricity defined by sjvxjvy. Equa¬ 
tion (I44l) . in this case, gives the same conductivities at 
fj, — En as in the isotropic case renormalized by the 
factor Vx/vy (vy/vx) for Gxx (^yy), which coincides with 
the conductivities obtained by the quasi-classical argu¬ 
ment. Note that in the case of the anisotropy induced 
by Vx ^ Vy, the ratio Uxxlv^yy = v^jv'^ is independent 
of the magnetic field (zero-field case was considered in 
Ref. 113), while if the anisotropy is produced by the tilt¬ 
ing, (Jxxivryy changes from sin“^ 77 /( 77 A) in zero field to 
A“^ in the high-held limit. The magnetic-Held-dependent 
anisotropy of the in-plane conductivity which saturates 
in the quantum limit can be used as an evidence of the 
tilting in Dirac electron systems with y close to the Dirac 
point. 
























X 

Figure 5. (Color online) Normalized conductivities CT|| = 
A“^(T|| and dx = Actx at /r = 0 with rj = 0.4 and 0.5 cal¬ 
culated from Eqs. and (HSI numerically using N = 150 
Landau levels, (a) The conductivities in the constant broad¬ 
ening approximation as functions of Tjhu}*; rj — 0 case is 
indicated by the dotted-dashed line, (b) The conductivities 
in SCBA as functions of xr with Ec/huj* = 100. In each fig¬ 
ure, the upper (lower) dashed line corresponds to the limiting 
value sin“^ v/i^v) zero field for rj — 0.5 {rj = 0.4).— 


IV. CONDUCTIVITY AT FINITE 
TEMPERATURES 

In the following, we discuss the magnetic-field depen¬ 
dence at finite temperature and the temperature depen¬ 
dence of the conductivity and the resistivity at /r = 0 
which can be directly compared with experiments 
The experimental measurements of the magnetic-field 
dependence of the in-plane conductivity in a-(BEDT- 
TTF) 2 l 3 organic conductor demonstrated a two-step 
decrease of the conductivity with increasing the mag¬ 
netic field— 44 In Ref. IH this two-step decrease was in¬ 
terpreted in terms of quasi-classical two-carrier model, 
where Dirac carriers and massive carriers coexist. Sim¬ 
ilar two-step behavior of the conductivity was also re¬ 
ported in Ref. [13. However, in Ref. [H detailed analysis 
of the magnetic-field dependence with fixed T revealed 
the existence of a novel minimum in the conductivity 
at a moderate magnetic field which scales as in the 
temperature range of 1.5 K < T < 5 K. In this sec¬ 
tion, taking account of the recent experiment;^ we show 



Figure 6. (Color online) Magnetic-field dependence of the 
conductivity at /r = 0 in the constant broadening approxima¬ 
tion for several temperatures without the Zeeman interaction 
(a) and with the Zeeman interaction included (b). 


an alternative explanation for the two-step behavior in 
terms of one carrier scenario of Dirac electrons. Fur¬ 
thermore, in the present explanation, the minimum of 
conductivity can be understood. At first, we consider 
the constant damping approximation which is frequently 
used in the literature— 44444444 Secondly, we show how 
the results obtained in the constant damping approxima¬ 
tion are modified if SCBA for the impurity scattering is 
used. We start with analyzing r] = 0 case in detail, and, 
after that, the effect of the tilting is considered. 

Hereafter, parameters are chosen with respect to the 
application to a-(BEDT-TTF)2l34 vp = 5 x 10^ m/s, 
r/A:^ = 3 K, and g = 2, if these values are not indicated 
explicitly. 


A. Conductivity in constant damping 
approximation 

In this section, we consider a simple case where the 
self-energy has the form E(e) = zF, with F being a phe¬ 
nomenological parameter. Magnetic-field dependence of 
the conductivity calculated from Eqs. JT]) and (fT^ for 
ji = 0 and 77 = 0 is shown in Fig. |6l Figure |6Ka) shows 
the result without the Zeeman interaction and IHKb) with 
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the Zeeman interaction. First, we discuss the case with¬ 
out the Zeeman interaction [Fig.IHlJa)]. We interpret that 
the decrease of a is associated with the crossover from the 
low magnetic-field region where Landau levels overlap to 
the quantum region where Landau levels are well sepa¬ 
rated. Actually, this crossover occurs when hujc ^ ksT. 
In the limit of B —>■ oo, the conductivity saturates at the 
value determined by the n = 0 Landau level. 

Next, we discuss the temperature dependence. The 
effect of the temperature on the conductivity is the op¬ 
posite in low and high magnetic fields. In the low-field 
region, increasing of the temperature enhances the con¬ 
ductivity. This will be due to the thermal activation 
of the carriers. In contrast, in the quantum region, in¬ 
creasing the temperature reduces the conductivity. This 
is because the n = 0 Landau level has a temperature 
smearing and, as a result, reduces the conductivity. 

At low temperatures, the crossover in Fig. [5] (a) be¬ 
tween the low magnetic-field region and high magnetic- 
field region can be discussed analytically from Eqs. © 
and m- We can use the Sommerfeld expansion for 
the Fermi-Dirac distribution —df/de = S{€ — fi) + 
(7r^fc^T^/6)d"(e — /x) which gives 

where cr"(0) = (Pcr{0)/d{fj,/T)'^ defined at /x = 0. cr"(0) 
can be found analytically from Eq. m 

cr"(0) = -2[l+ 47^ -p 87^ - 87V'(7^)] ■ (54) 


In large magnetic field (7 <C I), cr"(0) —>■ —2, and the 
conductivity in Eq. (1551) approaches ^{l—Tr'^kgT'^/ST^). 
This explains the increase of cr{T) as a function of T in 
high magnetic fields. When the magnetic field decreases, 
(t(T) reaches the universal value e^/(7r/i) at hjjJcjT k, 1.30 
which gives zero for the right-hand side in Eq. O, 
and roughly represents the crossover between low-field 
and high-held regions. At low magnetic helds (7 1), 

cr"(0) ~ I — and the behavior of the conductivity 

is given by 


a{T) 


tt/i [ 9F2 V 5 r4 / 


(55) 


This explains the increase of cr(T) as a function of T in 
low magnetic held. 

Next, we discuss the effect of the Zeeman interaction in 
Fig. El (b). When we include the Zeeman term there ap¬ 
pears the second step decrease in the high magnetic held 
due to the opening of the gap at /x = 0. The magnetic- 
held dependence of the conductivity in this region can be 
understood as follows. At zero temperature, the conduc¬ 
tivity per spin near /x = 0 in the spin-splitting region can 
be approximated as 





Figure 7. (Color online) Temperature dependence of the re¬ 
sistivity at /X = 0 for several magnetic fields. The dashed line 
shows the resistivity at zero magnetic field. The inset shows 
the position of the maximum fesTmax as a function of hwc 
with 3 = 0. 


At ^ = 0 


a(0,0) = 


tt/i 1 -I- {g^BB/2TY 


(57) 


This expression indicates that cr(0,0) decreases when 
ggsB > 2F. At hnite temperatures, cr{T) is obtained 
by Eq. 0. In this case, cr(T) starts to decrease for suf- 
hciently large magnetic held 2gBB > kBT. 

Finally, we discuss the temperature dependence of the 
resistivity. Figure [7] shows p(T, 0) for several magnetic 
helds. The characteristic feature of resistivity in Fig. [7] is 
the maximum in the temperature dependence which sep¬ 
arates the high-temperature region (/csT kuc) from 
the low-temperature region {kBT <C hujc)- In the high- 
temperature region, the resistivity is determined by the 
thermal activation of the carriers, which leads to the 
increase of the resistivity with lowering the tempera¬ 
ture. In the low-temperature region, where the n = 0 
level becomes isolated, the behavior of the resistivity is 
the opposite, since increasing the temperature reduces 
the density of states at /x = 0 due to the temperature 
broadening. Therefore, resistivity takes the maximum 
at fcaTinax ^ hujc- The result of the numerical cal¬ 
culation of Tmax on the {hiOc, kBT)-plaiie is presented 
in the inset in Fig. [71 which shows that at high fiwc, 
fcsTmax ~ 0.21{fujJc — 0.5r). In low magnetic fields, 
the maximum disappears when HlOc ~ 1.30 F, which also 
gives zero for the right-hand side in Eq. (P|) . 


B. Conductivity in SCBA 

In this section, we discuss how SCBA modifies the re¬ 
sults obtained in the previous section. The main fea¬ 
ture in SCBA is that the level broadening F(i3) increases 
as '/Bj^ As a result, the behavior of a changes when 
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Figure 8 . (Color online) (a) Magnetic-field dependence of 
the conductivity at /r = 0 in SCBA for several temperatures 
with >c = 0.02 and Ec = 0.05 eV. The inset shows vs 
temperature, (b) Magnetic-field dependence of fiojc, r(B), 
HbB, and the conductivity in the unit of e^/vr/i at T = 2.5 K 
(dashed line) in the same plot. The horizontal dotted-dashed 
line corresponds to T = 2.5 K. 


r(B) ^ ksT. This leads to a minimum of the conduc¬ 
tivity which is the central result in this section. 

FigurellKa) shows the conductivity for ^ = 0 and 77 = 0 
with the Zeeman interaction as a function of the mag¬ 
netic field in SCBA for several temperatures, where the 
self-energy was calculated numerically from Eq. (l30ll . In 
contrast to the simple two-step decrease in Fig.jSKb), the 
conductivity in SCBA takes a minimum in the moderate 
magnetic field region. We find that this minimum is due 
to the crossover between the “temperature-dominating” 
[r(i?) < ksT] and “impurity-dominating” [ksT < r(B)] 
regimes. In order to explain this, we show the energy 
scales of our problem, i. e., hwc, r(i?) [Eq. fisB, 

and ksT (dotted-dashed line) in Eig. [ 8 l[b). The dashed 
line is the conductivity for T = 2.5 K as a function of B. 
We can see that the minimum of a occurs at r(i3) ~ ksT. 
In the low-magnetic field region, T{B) is smaller than 
ksT, which corresponds to the temperature dominating 
region. In this case, a decreases according to Eq. (1551) . 
On the other hand, when T{B) > ksT, the conductivity 
slightly increases. This is because a can be approximated 
from Eqs. dial) and (l54ll as tr = l—Tr‘^kgT'^/3r'^{B) which 
approaches 1 with increasing r(i?). Then, in the high- 



Figure 9. (Color online) Magnetic-field dependence of the 
conductivity at 77 = 0 in SCBA for T = 4.5 K, x = 0.02, 
and Ec = 0.05 eV calculated from Eqs. dm and (130 II (solid 
line). Closed squares (triangles) show a± (cr||) with r; = 0.5 
calculated from Eqs. (0 and da [Eq. dmi using summation 
over 30 Landau levels. The horizontal axis is renormalized as 


magnetic field region, T{B) starts to decrease, and, as a 
result, (T decreases due to the Zeeman interaction. 

In order to confirm this interpretation, we calculate 
Bmin of various temperatures numerically. As shown in 
the inset of Fig. UKa), i?min scales as in agreement 
with the relation T{B) ^ \/B ^ ksT. The minimum 
of a disappears at high temperatures [e. g., “T = 9 K” 
curve in Fig.[5l[a)]. This is because ksT > T{B) holds in 
all range of magnetic fields. 

In the high-magnetic-field region where the Zeeman 
term reduces a, the B dependence of a is different from 
that obtained in the previous section. At T = 0, accord¬ 
ing to Eq. (l38ll . we obtain 


a(0,0) 




for B < Be, 

fo, 


otherwise. 


(58) 


where Be = >cehvp/{g^fM^). At finite temperatures, the 
conductivity in the region B > Be decreases exponen¬ 
tially with magnetic field, cr{T) ^ exp(—/rsB/ZcsT). 

Next, we study how our results for the magnetic-field 
dependence of the conductivity are modified by the tilt¬ 
ing of the Dirac cone. First of all, there are 77 -dependent 
prefactors of a± and a\\ [Eqs. ([33] and (HUl) ]. We will 
show that the magnetic-field dependencies of (t_l and 
are not affected so much by 77 except for these pref¬ 
actors. Figure [3] shows the normalized conductivities 
a± = Xa± and (T|| = A“^(T|| as a function of B for 
T = 4.5 K and x = 0.02 where there is no Zeeman 
term. Note that the magnetic field is also renormalized 
by 77 as huje X~^^^fujJe- Thus, the horizontal axis in 
Fig. Ois .8 = B{1 — 77 ^)“^/"^. Closed squares (triangles) 
in Fig.inishow a± (o-||) calculated from Eqs. ([T]) and (fT3ll 
[Eq. ([TKll ] using numerical summation. The solid line is 
the same curve as in Fig. [5][a) for T = 4.5 K which is 
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Figure 10. (Color online) Resistivity at /i = 0 as a function of 
temperature for several magnetic fields calculated in SCBA 
with jr = 0.02 and Ec = 0.05 eV. 


shown for comparison in the case of 77 = 0. Figure [5] 
indicates that 77 has very small effects on the normal¬ 
ized conductivities. Only in the weak magnetic-field re¬ 
gion, CTii is slightly larger than d'±. This difference can 
be understood in terms of Eqs. dMl) and H40p . In this 
case, there are Landau level mixings, and the contribu¬ 
tion from higher Landau levels for (Tj_ is smaller by the 
factor than that for i 7 ||. 

The physical reason of this small rj dependence will 
be as follows. The basic elements of our explanation for 
the magnetic-field dependence are the Landau level quan¬ 
tization and the magnetic-field dependence of the level 
broadening. This explanation remains qualitatively the 
same even if the tilting is introduced. Moreover, in the 
quantum limit, the conductivity is solely determined by 
u = 0 Landau level. In this case, n = 0 Landau level is in¬ 
sensitive to 77 except for the rescaling factors A and A“^. 
This will be the reason why the tilting has very small 
effect on the magnetic-field dependence of the conduc¬ 
tivity. Note that in the high-magnetic field region, where 
o’-L (cii) reduces due to the Zeeman splitting of n = 0 
Landau level, the 77 dependence comes only through the 
prefactor A“^ (A). 

Finally, we briefly mention the temperature depen¬ 
dence of the resistivity in SCBA shown in Fig. [TU] for 
several magnetic fields. The explanation of this behavior 
remains qualitatively the same as in the previous section 
for the constant broadening approximation. The main 
difference of the present case from Fig. [7] is that the max¬ 
imum in the temperature dependence of the resistivity 
survives at small magnetic fields and becomes more pro¬ 
nounced as magnetic field decreases. This behavior can 
be understood from Fig.[7]which indicates that the n = 0 
level becomes sharper when magnetic field goes to zero, 
since T{B) for >c = 0.02 becomes exponentially small at 
B = 0. 


V. SUMMARY 

In summary, we have made an analytical calculation of 
the longitudinal conductivity of two-dimensional mass¬ 
less fermions with a tilted Dirac cone in the framework 
of the Kubo-Basting-Stfeda linear-response formula in 
strong magnetic field. The main difference of the present 
case from the case of the isotropic Dirac cone is that 
the tilting introduces Landau level mixing which breaks 
usual selection rules for the matrix elements of velocity 
operators which leads to the existence of matrix elements 
between the Landau levels with arbitrary n. 

First, we have analyzed the conductivity at zero tem¬ 
perature in the quantized magnetic field limit. In this 
case, we found that the conductivities at ^ are 

given by axx = + 2 Y — Nr]'^) and Uyy = 

+ 2N). The factors A and A“^ are purely geo¬ 
metrical, due to the ellipticity of the quasi-classical elec¬ 
tron orbit. Apart from these factors, ayy remains the 
same as in the case without the tilting, while axx is re¬ 
duced by Nr]^. We have explained this result as fol¬ 
lows. The conductivity at /i = En is determined by the 
quantum-mechanical average of (Ar;)^ calculated with 
respect to the iVth Landau level wave function. If the 
tilting of the cone is in the ky direction in the energy- 
momentum space, the center of the Nth eigenfunction is 
displaced hy xn along the x direction in the real space. 
We have proven that this displacement gives the reduc¬ 
tion of (Aa:)^, while (Aj/)^ remains unchanged. We note 
that this fact is purely quantum effect and has no quasi- 
classical interpretation. 

Second, we studied how the zero-temperature conduc¬ 
tivity at /J, = 0 depends on magnetic field and impurity 
scattering. Using numerical calculations, we have demon¬ 
strated that axx at /i = 0 remains independent of mag¬ 
netic field and short-ranged disorder, which is similar to 

the case without tilting, while ayy shows a crossover from 
2 • — 1 

the value ^ in zero magnetic field to the value 

if in the strong-field limit. We propose that this 
fact can be used as an experimental evidence of the tilting 
of the Dirac cone. Indeed, there are two possible sources 
of the anisotropy of the conductivity, namely, owing to 
the difference in the Fermi velocities Vx ^ Vy, and due to 
the tilting. In the former case, ayy j axx = v^lvy does not 
depend on magnetic field, while in the latter case ayy j axx 
changes from yT^^^sin”^ 77/77 in zero magnetic field to 
l — if' in the high-held limit. For a-(BEDT-TTF) 2 l 3 with 
77 Ri 0.8, ayy/axx changes, respectively, from 0.69 to 0.36 
which can be detected experimentally. The magnetic- 
held dependent anisotropy of the in-plane conductivity 
can be regarded as a characteristic feature of a tilted 
Dirac cone. 

Next, we have studied the conductivity at hnite tem¬ 
peratures where we set 77 = 0 at the beginning. We have 
demonstrated that as a function of the magnetic held 
the conductivity at the Dirac point undergoes a two-step 
decrease where the hrst decrease is associated with the 
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crossover to the quantum limit and the second decrease is 
due to the spin splitting of the n = 0 level. For the effect 
of the impurities, the results are given for two models: 
the constant broadening approximation and SCBA. The 
analysis based on the SCBA shows that the broadening 
of the Landau levels increases with the magnetic field as 
\/B, which leads to the minimum in the magnetocon¬ 
ductivity at ksT ^ T{B) ^ which separates the 

temperature dominated region ksT > r(B) from the im¬ 
purity dominated region ksT < r(i3). We have demon¬ 
strated that this behavior remains qualitatively the same 
even in the presence of the tilting of the cone. 

The results of Sec m can be directly compared 
with recent experiments in q:-(BEDT-TTF) 2 I 3 organic 
conductori^i^ We note, that in contrast with Ref. 1^ . 
our analysis is able to explain the magnetic-field depen¬ 
dence of the conductivity in a-(BEDT-TTF) 2 l 3 using 
only one type of Dirac carriers. Moreover, the minimum 
in the magnetoconductivity, which scales with tempera¬ 
ture as , can be regarded as a strong evidence in favor 
of the one-carrier scenario of the magnetotransport. 

In conclusion, we comment on the minor importance of 
the tilting of the Dirac cone for understanding the exis¬ 
tent experiments in magnetic field in q;-(BEDT-TTF) 2 I 3 
near /i = 0. We note that the key ingredients of our 
explanation of the magnetic-field dependence of the con¬ 
ductivity, namely, the Landau quantization, the existence 
of the zero mode, and the magnetic-field dependence of 
the level broadening remain qualitatively the same even 
in the presence of the tilting. That is the reason why 
our main results in Sec. El such as two-step behavior of 
the conductivity in magnetic field and the minimum in 
the magneto-conductivity remain valid in the presence 
of the finite tilting. In order to detect the effect of the 
tilting on the in-plane transport, the authors propose to 
measure the magnetic-field dependence of Oxxj^yy This 
effect survives at finite temperatures and can be used as 
a marker of a tilted cone. 
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Appendix A: Matrix elements of velocity operators 


The matrix elements of the velocity operators in 
Eqs. (fTTIl and (fT^ for the Landau levels in Eq. Q can 
be expressed in the following form: 

{m\vx\n) = iXCnm - sgn{n)F^^^'> , (Al) 

{m\vy\n) = Cnm + sgn{n)F^]^'> + sgn{m)F^°^'> 

- r]sgn{m)sgn{n)F^]^'> , (A2) 

where Cnm = (2 - dno)~^^^{2 - 

/ OO 

-OO 

i^j = 0,1, and \n\ > i, |m| > j. The coefficients Fnm 
can be calculated using the formulap^ 

/ OO 

Hm{x + y)Hn{x + z)dx 

-OO 

= 2"7ri/2m!z"-™L;)r'"(-22/z), {m < n), (A4) 

which gives 

\/R 


p(10) ^ ATm Vl^l A2 




Jml 


p(01) _ ATmV^r^j-\n\-\m\ + lf 2 a2 \ 
^nm n n -r — l ^ 

^J\m\ ' ' 


(A5) 

(A6) 

(A7) 


f(11) = a2 n 

^ nm y \ \-^\m\ — l \ l ‘-^n,mh 


(AS) 


where 


'™'exp^-yA2^^^ , (A9) 


and we used the property of Laguerre polynomials with 
integer n and m: 


{-xY 


■LI 


(AlO) 


Substituting Eqs. (IA5])-(IA81) in the matrix elements (ED 
and (IA2I) and using the recurrence relations for Laguerre 
polynomials 

- CFFCh (All) 


dm| 


|m| —1 




dm| 


|m| 


we obtain the expressions (ED and ED- 


(A12) 
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Appendix B: Conductivity in the 77 — >■ 0 limit 

If 77 = 0, the matrix elements in Eqs. (ED and m become 

{m\vx\n) = -iCnm [sgn(n)(5|„||m|+i - sgn(m)(5|m||„|+i] , 
{m\vy\n) = Cnm [sgn(n)(5|„||m|+i + sgn(m)i5|m||„|+i] , 


(Bl) 

(B2) 


where Cnm = (2 — dno) ^^^(2 — i5mo) In this case, Eqs. (Ill)) and (IT5|) for the conductivities reduce to Eq. (fT71) . 
In order to make a summation over n in Eq. (ED, we make a summation over a and a' with the help of the identity 


(x — a^/n)^ + 7 ^ 2i \x — a^/n — 77 x — a^/n + 77 


(B3) 


which yields 


o-Q 




7+ , T'- 


4^\n+l + 7^ 77 + 1 + 77/ Vn + 7+ 77 + 71 


7+ , 'I'- 


(B4) 


where we introduced a shorthand notation 7 ± = (7 ± ix)/huJc- Now the summation over n can be easily done with 
the help of the formula 


E 

n=0 


1 


1 


(77 + a) (77 + 6 ) b — a 


[^{b) - - 0 ( 0 )], 


(B5) 


where ip denotes digamma function. The final result is 


Mx) = 2 


7^+7^+ ( 7 ^-7^)^ _ 7+7- ( 7 ^ - 7^) 7.2 ) _ ^7 2 

+ 27+7-(1-(7^-7^)^) l-( 7 i- 7 ^)^ 


(B 6 ) 


which gives Eq. (ED- 

In the limit ijAn^m —>■ 0 we can expand Eqs. (fTTll and (TT^ in powers of 77 . Up to the order of 77 ^, Vx and Vy have 
matrix elements up to | 777 | = | 77 | ± 3. The expansion has the form 


(77l|Va;|77) jXCnm ^ ^ ^|m| = |n|+j ^\n\ = \m\+j 


o{v% 




(B7) 


where 


01 ) 


02 ) 


03) 


= Sgn(777) j 1 “ 2 + 1 “ Sgn(77777)A/E^ (|77| + 1)^^^ - Sgn(77777) |77p/^ 


= y/\n\ + 1 ( V|n| + 2 - sgn(77777)-\/E) 77 , 


= 1 Sgn(777)^/ (|77| + 1)(|77| + 2) ( V| 77 | + 3 - Sgn(77777) , 


(B 8 ) 

(B9) 

(BIO) 


and 


(rTT-lUylTT.) — A Cnm^^ ^'^|m| = |n|+j + '^|ra| = |m|+i +10(7 ), 

3=0 


(Bll) 


where 


&()’) = - 27777 , 


bn^ = Sgn(m) |l - ^ (V|77| + 1 - Sgn(77777) 


3 r 


(|77| + 1)^^^ + Sgn(77777)|77p/^ 


(B12) 

(B13) 

(B14) 
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Appendix C: Calculation of ((Ar;)^)/ 


In the Landau gauge, which is used in this paper, the x operator is bounded, and calculation of the following matrix 
elements is straightforward. 

{^kN\ (x - Xjv)^ l^'fe'Ar) = 2-Sj^ 1/ dxhf^^{x)xP 

/ OO nCO >, 

dxhiN\{x)hiN\_i{x)xP + {1 - Sm) dxh1j^^_^{x)xP > (Cl) 

-OO j —OO J 

where p = 1,2, and the integrals are calculated as 

/ OO r nOO 

dxh^N\{x)xP = ^ {2\N\ + 1), / dxhiN^{x)hiN\-iix)xP = 6. 

-OO J — OO 



(C2) 


which justifies Eqs. (PI and (1481) . However, the y operator is unbounded, which usually means that the y integral 
is ill defined and depends on the choice of the limits. We choose the following regularization. We add exp 
where e —>■ 0+, and, after that, expand the integration limits to infinity: 


poo 

r , 1 

exp 

{k—k')^ 

e2 

/ dy exp 

y —OO 

i{k-k)y- —y^ 

= 27r 



dy exp 

X) 

dy exp 


i{k-k')y- —y^ 


i{k- k')y- —y" 


y = —2m 


d 


{k—k'y 

- 


y^ = -2 tt 


dk 

^2 exp - 


dA:2 


\/^e 


2’Kdik — k'), 

—2m6'{k — k'), 
-2Tr6"{k-k'). 


(C3) 

(C4) 

(C5) 


The matrix elements of the y-operator correspond to the integrals of total derivatives of h'^^^{x) and h\N\{x)h\pf\_i{x), 
which give zero. The matrix elements of y^ can be expressed as 


{'i’kNly'^l'^k'N) = 


2 — Sno 


dx 


2sgn(n)?7 / dxh\f^^{x)h\j^^_-y{x) + (1 - 5no) [ dx (h\ff^_^{x)] 1 , (C6) 

J —OO J —OO ' ' ) 


where 


/ \ 2 ^ 

dx[h\j^^{x)j =-(2|iV| + l), and J dxh\j^^{x)h\j^^_^{x) = 0, 


(C7) 


which yields Eq. (I50|. 


* iprosk@hosi.phys.s.u-1okyo.ac.jp 

^ H. Fukuyama, Y. Fuseya, M. Ogata, 
A. Kobayashi, and Y. Suzumura, 

Physica B: Condensed Matter 407, 1943 (2012) 

^ K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, 
M. I. Katsnelson, I. V. Grigorieva, S. V. Dubonos, and 
A. A. Firsov, Nature 438, 197 (2005) 

® X.-L. Qi and S.-C. Zhang, 

Rev. Mod. Phys. 83, 1057 (2011) 

^ N. Tajima, S. Sugawara, M. Tamura, Y. Nishio, and 
K. Kajita, J. Phys. Soc. Jpn. 75, 051010 (2006) 


® S. Katayama, A. Kobayashi, and Y. Suzumura, 
J. Phys. Soc. Jpn. 75, 054705 (2006) 

® H. Kino and T. Miyazaki, 

J. Phys. Soc. Jpn. 75, 034704 (2006) 

^ A. Kobayashi, S. Katayama, Y. Suzumura, and 
H. Fukuyama, J. Phys. Soc. Jpn. 76, 034711 (2007) 

® K. Kajita, Y. Nishio, N. Tajima, Y. Suzumura, and 
A. Kobayashi, J. Phys. Soc. Jpn. 83, 072002 (2014) 

® P. Richard, K. Nakayama, T. Sato, M. Neupane, Y.- 
M. Xu, J. H. Bowen, G. F. Chen, J. L. Luo, N. L. 
Wang, X. Dai, Z. Fang, H. Ding, and T. Takahashi, 
























15 


Phys. Rev. Lett. 104, 137001 (2010) 

“ T. Osada, J. Phys. Soc. Jpn. 77, 084711 (2008) 

T. Morinari and T. Tohyama, 

J. Phys. Soc. Jpn. 79, 044708 (2010) 

I. A. Luk’yanchuk, A. A. Varlamov, and A. V. Kavokin, 
Phys. Rev. Lett. 107, 016601 (2011) 

I. Proskurin and M. Ogata, 

J. Phys. Soc. Jpn. 82, 063712 (2013) 

T. Konoike, M. Sato, K. Uchida, and T. Osada, 

J. Phys. Soc. Jpn. 82, 073601 (2013) 

N. Tajima, S. Sugawara, R. Kato, Y. Nishio, and K. Ka- 
jita, Phys. Rev. Lett. 102, 176403 (2009) 

N. Tajima, M. Sato, S. Sugawara, R. Kato, Y. Nishio, and 

K. Kajita, Phys. Rev. B 82, 121420 (2010) 

S. Sugawara, M. Tamura, N. Tajima, R. Kato, 

M. Sato, Y. Nishio, and K. Kajita, 

J. Phys. Soc. Jpn. 79, 113704 (2010) 

N. Tajima, T. Yamauchi, T. Yamaguchi, M. Suda, 

Y. Kawasugi, H. M. Yamamoto, R. Kato, Y. Nishio, and 

K. Kajita, Phys. Rev. B 88, 075315 (2013) 

M. O. Goerbig, J.-N. Fuchs, G. Montambaux, and 

F. Piechon, Phys. Rev. B 78, 045415 (2008) 

M. O. Goerbig, J.-N. Fuchs, G. Montambaux, and 

F. Piechon, Europhys. Lett. 85, 57005 (2009) 

T. Kawarabayashi, Y. Hatsugai, T. Morimoto, and 

H. Aoki, Phys. Rev. B 83, 153414 (2011) 

J. Sari, G. Toke, and M. O. Goerbig, 

Phys. Rev. B 90, 155446 (2014) 

M. Trescher, B. Sbierski, P. W. Brouwer, and E. J. 
Bergholtz, Phys. Rev. B 91, 115135 (2015) 

V. Lukose, R. Shankar, and G. Baskaran, 
Phys. Rev. Lett. 98, 116802 (2007) 


N. H. Shon and T. Ando, 

J. Phys. Soc. Jpn. 67, 2421 (1998) 

N. M. R. Peres, F. Guinea, and A. H. Castro Neto, 

Phys. Rev. B 73, 125411 (2006) 

S. G. Sharapov, V. P. Gusynin, and H. Beck, 

Phys. Rev. B 67, 144509 (2003) 

V. P. Gusynin and S. G. Sharapov, 

Phys. Rev. B 71, 125124 (2005) 

T. Morinari and T. Tohyama, 

Phys. Rev. B 82, 165117 (2010) 

Y. Suzumura, 1. Proskurin, and M. Ogata, 

J. Phys. Soc. Jpn. 83, 023701 (2014) 

Y. Suzumura, 1. Proskurin, and M. Ogata, 

J. Phys. Soc. Jpn. 83, 094705 (2014) 

A. Kobayashi, Y. Suzumura, and H. Fukuyama, 

J. Phys. Soc. Jpn. 77, 064718 (2008) 

M. Monteverde, M. O. Goerbig, P. Auban-Senzier, 

F. Navarin, H. Henck, C. R. Pasquier, C. Meziere, and 
P. Batail, Phys. Rev. B 87, 245110 (2013) 

T. Ozawa, T. Yamauchi, N. Tajima, 

R. Kato, Y. Nishio, and K. Kajita, 

Meeting abstracts of the Physical Society of Japan 69, 862 (2014) 
T. Morinari, T. Himura, and T. Tohyama, 

J. Phys. Soc. Jpn. 78, 023704 (2009) 

A. Bastin, C. Lewiner, O. Betbeder-Matibet, and 
P. Nozieres, J. Phys. Chem. Solids 32, 1811 (1971) 

P. Stfeda and L. Smrcka, 

Phys. Status Solid! B 70, 537 (1975) 

S. G. Sharapov, V. P. Gusynin, and H. Beck, 

Phys. Rev. B 69, 075104 (2004) 

1. S. Gradshteyn and I. M. Ryzhik, Table of integrals, se¬ 
ries, and products, 7th ed. (Elsevier/Academic Press, Am¬ 
sterdam, 2007). 


25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 



